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1.  INTRODUCTION 


The  solution  of  the  wave  equation  constrained  to  a 
cylindrical  geometry  gives  rise  to  a  radial  dependence  in  terms 
of  Fourier-Bes sel  functions.  This  analysis  is  that  appropriate  to 
an  annular  duct  and  therefore  the  eigenfunctions  include  both 
Bessel  and  Neumann  contributions.  At  any  radius,  the  argument  of 
the  Fourier-Bessel  function  is  determined  by  an  eigenvalue  which 
depends  on  both  the  radial  node  number  (1)  and  the  order  of  the 
eigenfunction  (n).  These  eigenvalues  are  determined  from  the 
boundary  condition  of  zero  derivative  at  the  inner  and  outer 
radii  of  the  duct. 

The  method  used  herein  seeks  to  construct  solutions  of  the 
large  order  Fourier-Bessel  functions  from  those  of  the  zeroth 
order  by  obtaining  coefficients  to  a  truncated  infinite  series 
expansion.  The  coefficients  are  eigenvectors  which  are  associated 
with  modified  versions  of  the  large  order  eigenvalues. 

The  number  of  terms  necessary  for  accuracy  of  the  truncated 
infinite  series  is  investigated  and  the  resulting  large  order 
eigenfunctions  are  plotted. 


2 . THEORY 


2.1  Basic  Equations 


The  solution  of  the  wave  equation  (with  or  without  mean 

f  low) 


[  1  i_  l¥(r,t)  =  0 

'  a*  at*  ' 

where  flr^t)  is  a  wave  function  and  a.  is  the  speed  of  the  wave 
motion,  results  in  a  radial  dependence  given  by 

d*R?  +  1  dR?  +  fkT"  -  n*  '|RT  =0 
dr'  r  dr  k  r* ' 

where  RJ ( r )  are  a  complete  set  of  orthonormal  radial 
eigenfunctions  called  Fourier-Bessel  functions  and  RJ(r)  e  R„ ( k£ 
r),  i.e.  n  is  the  order  of  the  function  and  1  is  the  radial  node 
number  ( 1=0 , 1 , 2 , . . . ) .  Specifically,  the  RJ(r)  are  a  combination 
of  Bessel  and  Neumann  functions.  The  Fourier-Bessel  function 
differs  from  the  Bessel  function  because  it  depends  on  n  which  is 
constrained  by  the  number  of  circumferential  wavelengths  fitted 
around  the  duct,  i.e.  n  is  the  circumferential  mode  number. 

Hence,  because  the  circumferential  dependence  is  harmonic,  the 
"Bessel"  function  has  a  Fourier  component  and  therefore  we  call 
it  a  Fourier-Bessel  function.  Solutions  of  this  equation  are 
admitted  only  for  certain  values  of  the  eigenvalue  kj  which  is 
determined  from  the  boundary  condition 
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dRJ  =  0 
dr 

at  r=1.0,h  where  h  is  the  ratio  of  the  huh  of  the  duct  to  the 
outer  radius.  This  boundary  condition  arises  because  the  radial 
gradient  of  the  wave  function  is  zero  at  the  inner  and  outer  duct 
radii;  for  example,  the  radial  velocity  must  vanish  at  the  hub 
and  the  tip  for  non-porous  walls. 

Normalisation  of  the  eigenfunctions  fixes  their  magnitudes 
absolutely  and  is  achieved  by 


f1  rR£  <r)l£  (r)  dr  =  6^ 

2.2  Expansion  Of  The  Higher  Order  Modes 

Following  the  approach  of  Namba  ( ref . 1 )  we  expand  the  nth 
order  eigenfunction  in  terms  of  the  zeroth  mode 

R?  ( r )  =  *E  BTR*  ( r  > 

*  m>m  0  * 

where  L  is  as  large  as  we  wish.  Hence,  we  have 

r*‘  +  1  b£.  r:'  +  (kf‘  -  nM  BJ,  R*  *  0 

"*  r  v  r»  y 

and  also 

+  1  Ri'  +  KK,  =  0 

r 

where  summation  over  m  is  implied  for  both  equations. 

-  i  -  js)b;.h;  =  o  . 

V  ?  n* ' 

Putting  KJ  =  In  ,  then  multiplying  by  rRj (r)  and  integrating 
|*  (KJ*  -  £  J^Rjrdr  -  £  X  iCR/dr  =  0 

Bi  S  =  (  R^j  +  ift«wj  )BL 

where  ^  =  '[  1  Ri(r)Rt  (r)  dr  (2.2.1) 

t!  r 

and  KJ  and  are  respectively  the  eigenvalues  and  eigenvectors 
of  a  matrix  " 
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(2.2.2) 


Rj  +  iSj  . 

n 

Note  that  j'  rRJ  RJ  dr  =  j‘  B£.  R^R^r  dr 

<v  =  ^  S-U 

£  ^  =  1  (2.2.3) 


i.e.  the  normalisation  condition  applied  to  the  eigenfunctions 
means  that  the  sum  of  the  squares  of  the  eigenvectors  is  one. 


2.3  The  Boundary  Condition 

Since  R£ ( r )  is  a  combination  of  Bessel  and  Neumann 
functions,  we  can  represent  it  by 


R£ ( r )  =  AJ£  <  r )  +  BY* ( r ) 

where  A  and  B  are  constants  depending  on  n  and  1.  We  have  dR J =  0 
at  r=1.0,h  dr 

hence,  k£Aj;  (kjh)  +  kJBY^(k^h)  =  0 

and  k£Aj;<k^)  +  kj  BY«  ( kj  )  =  0  . 

Assuming  kj  *  0,  this  gives 


B  =  -  i;  ( k?  )  =  -  JA  (kih) 


(2.3.1) 


y;o<;) 


YA  ( kjh) 


Since  our  solution  does  not  extend  to  r=0,  B*0  and  therefore 


JA<k;  )HJ  (kjh)  -  JA  (kJh)Y;  (k^  )  =  0 


Abramowitz  (ref. 2)  gives  the  solution  of  this  as 


lC  h  =  B  +  E  +  q  -  p~  +  r  -  4pq  ■»  2p  +  ...  (2.3.2) 

p  "  e* 

where,  if  p  =  4n*  and  X  =  1.0 /h 

0  =  It  p  =  u  +  3 

X  -  1.0  8X 
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B  =  -  J.  (kf  ) 
A  Y,  ( kj  ) 


Hence,  for  the  zeroth  order  eigenfunction,  both  coefficients  can 
easily  be  determined  and  the  problem  now  becomes  one  of 
evaluating  the  integral  (2.2.2)  and  extracting  the  eigenvectors 
from  the  resulting  matrix. 


2.5  Zeroth  Node  Number  Terms 

R£(r)  is  a  complete  set  of  orthonormal  functions  and 
therefore  the  1=0  term  must  be  included.  Clearly,  the  boundary 
condition  formula  (2.3.2)  cannot  be  used  here  since  1=0  results 
in  a  singularity. 

Actually,  we  see  that  ki  =0,  because  n=l=0  means  that  there 
are  no  nodes  in  the  radial  and  circumferential  directions,  i.e. 
the  wavelengths  in  each  case  are  infinite,  and  therefore  the 
formula  cannot  be  derived.  If  kj=0,  the  arguments  of  the  Bessel 
and  Neumann  functions  are  zero.  However,  this  results  in 
Y*(0)  ->  -«*>  and  ;.  B  =  0  f or  k£  =  0.  Hence, 


r;  ( o )  =  aj;  (  0  ) 

for  n=0  R*  ( 0 )  =  A  since  Jl(0)  =  1.0 


i.e.  Rl  =  constant  . 


We  have  £  rRJdr  =  1.0 


S‘(l-h*  )  =1.0 
2 


Equation  (2.2.1)  gives 
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R^,  =  -  2 log,  h 
1.0-h*- 


We  are  now  able  to  form  all  values  of  the  integrand  of  R^  for 

Sr  J  =  0,1,2/  »  •  *  r  I/-  1  • 


2.6  Outline  Of  The  Method  Of  Solution 


Use  of  a^  Simpson's  rule  approximation  will  now  be  sufficient 
to  evaluate  for  all  values  of  m  and  j.  Addition  of  IcT/n*-  to 
the  diagonal  elements  of  the  matrix  thus  obtained  will  provide  us 
with  a  matrix,  the  eigenvalues  and  eigenvectors  of  which  will  be 
KJ*  and  Bj^  respectively. 

We  note  in  passing  that  the  matrix  must  be  real  and 
symmetric  and  therefore  its  eigenvalues  will  be  real  which 
implies  that  its  eigenvectors  are  real.  Therefore,  the  nth  order 
eigenfunctions  must  be  real.  Furthermore,  as  n  -»  the  influence 
of  the  k£/n*  becomes  less  and  less  and  hence  the  matrix  becomes 
constant  for  given  L  and  h,  i.e. 


RJ  ~  Rj1  +  0(l/nl)  . 


3.  RESULTS 


A  computer  program  based  on  this  theory  has  been  written 
which  demonstrates  the  behaviour  of  the  large  order 
eigenfunctions . 

It  is  useful  first  to  observe  from  fig.l,  which  is  taken 
from  ref.l,  that  the  eigenvalues  KJ*  of  the  matrix  given  by 
(2.2.2)  decrease  monotonical ly  as  n  increases,  but  that  the 
limits  are  all  greater  than  unity  and  differ  with  1. 

Additionally,  for  given  1,  KJ  can  be  made  closer  to  unity  by 
increasing  the  number  of  retained  terms  L.  This  behaviour  is  as 
expected  since  as  n-*»,  the  matrix  tends  to  TUj  ,i.e.  the  n 
dependence  drops  off  as  1.0/n *"  and  therefore  a  limit  is  reached. 
This  limit  is  not  unity  because  an  infinite  number  of  terms  are 
not  retained.  Namba  (ref.l)  and  McCune  (ref, 4)  discuss  this 
behaviour  in  more  detail. 

The  eigenvectors  E£,  associated  with  the  KJ*-  are  calculated 
subject  to  the  normalisation  condition  (2.2.3)  and  these  together 
with  the  zeroth  order  Fourier-Bessel  functions  provide  the  n*" 
order  Fourier-Bessel  functions  which  are  plotted  in  figs.  2,  3, 
and  4.  It  is  noticed  that  for  given  n  at  1=0,  most  of  the  effect 
of  the  radial  variation  is  concentrated  near  to  the  outer  radius, 
but  that  as  1  increases  the  peak  shifts  to  the  inner  radius.  If 
we  examine  the  eigenvalues  for  each  RJ  and  take  the  product  r 
then  we  can  choose  for  the  peaks  to  occur  at  Kir* =1.0  (where  rt 
is  the  radius  at  the  peak)  by  judicious  choice  of  L.  Hence,  as  1 
increases,  KJ  increases  and  therefore  the  peak  occurs  at  smaller 
r,  i.e.  towards  the  hub.  Conversely,  as  n  increases  at  constant 
1,  the  peak  will  move  to  the  outer  wall. 

The  reason  why  the  optimum  number  of  terms  for  the  expansion 
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(L)  was  chosen  such  that  KJ  =1.0  is  explained  in  the  next 
section.  It  should  be  noted  that  an  asymptotic  expansion  does  not 
necessarily  become  more  accurate  as  more  terms  are  taken.  The 
value  of  L  which  made  K£  cf  =1.0  turned  out  to  be  6  or  7. 


4.  DISCUSSION  OF  RESULTS 


Abramowitz  (ref. 2)  gives,  for  large  order  Bessel  and  Neumann 
functions 


J„  (nz )  -  (  4C  )*  Ai(n *C!  +  ... 

vl-z‘;  n*‘ 


Y„  ( nz ) - f  4  C,  Bi(n*C)  +  .  .  . 

^  1  ~  z**  '  n*» 


where  Ai  and  Bi  are  the  Airy  functions  given  by 
Aifif^)  ~  i  (n**4  sin  (^+H)  + 

Bi<n*£)  ~  1  (n^Ci^cos  fC+i'i  +  ... 

it*  4' 

for  C,< 0  and 

2/3(-Q)*  =  (z*■-l^-  cos“(l/z)  for  z>l. 

If  we  put  z=Kjr,  then  nz  becomes  nK£r  =  kj  r ,  i.e.  the  argument  of 
our  Fourier -Bessel  functions.  We  can  now  construct  an  expression 
for  the  large  order  Fourier-Bessel  functions:- 

R„(k?r)  ~  A  (  4C  (rftc  )~Vsin(C+Tt/4  )  +  B  cos  ( £+ir/4 )  1 
nv  'T-?"'  ir'h  '  A  ' 


or 

R„(kjr)  ~  A  (4  f*  (sin((,+*/4 )  +  B  cos(£+ir/4)) 
(nir)v*-'l-z*“'  '  A 


To  get  quantitative  results  from  this  formula  we  need  to 
find  KJ .  This  cannot  be  found  from  the  boundary  condition  (2.3.2) 
despite  the  fact  that  K£=  k£/n  because  this  formula  does  not  hold 
for  large  n. 

Furthermore,  the  expression  for  £  is  only  valid  for  £  real, 
i.e.  z>l.  If  this  is  not  the  case,  an  alternative  expression  for 
C  is  used  which  leads  to  exponential  terms  for  the  Airy  functions 
and  their  derivatives.  Interested  readers  are  referred  to 
Abramowitz  (ref. 2).  Here,  we  simply  observe  from  (4.1)  that 
K£r=l,  i.e.  z=l.  implies  a  peak  for  R£(r). 
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5.  CONCLUSION 


r  ^  It  is  clear  from  the  partial  derivation  of  a  formula  (4.1)' 

which  attempts  to  evaluate  analytically  the  large  order 
Fourier-Bessel  functions  that  problems  are  posed  such  that 
results  can  only  be  extracted  with  great  difficulty.  A  major 
hurdle  is  the  evaluation  of  the  eigenvalues  at  large  n. 

By  contrast,  an  infinite  series  expansion  truncated  after 
an  appropriately  large  number  of  terms  provides  the  behaviour 
hinted  at  (but  not  described)  in  the  formula  (-4-.  and  utilises  a 
simpler  approach.  The  evaluation  of  the  integral  (-2.2.2f  is 
central  to  the  method  and  this  presents  no  difficulty. 

Some  correlation  of  the  results  obtained  is  provided  by 
reference  to  the  peak  value  of  R£(r)  and  *ts  behaviour  with 
respect  to  changing  n  and  1.  From  the  examples  taken  6  or  7 
terms  are  seen  to  be  necessary  for  the  expansion.  Values  of  RJ(r) 
for  large  n  show  a  decreasing  dependence  on  this  parameter  and  in 
fact,  for  n-x»,  R?  ( r )  -*  Rj’(r)  which  depends  simply  on  R«:  given  by 
(2.2.1). 

^This  work  has  examined  the  behaviour  of  a  perhaps  recondite 
function  at  large  orders,  but  it  is  one  which  is  essential  for 
correct  analysis  of  the  radial  dependence  of  a  wave  function 
where  tjie  inner  boundary  of  a  cylindrical  geometry  has  a  finite 
value.  ;  f- («•.-'  •'  ,  ' .  -,r  •  f , , .  •  ;  .  .. 


S .  ROBERTS  ( S . 0 . ) 
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